27/01/95 
February 1, 2008 



Further solutions of critical ABF RSOS models 



Yu-kui Zhou[ 

Mathematics Department, University of Melbourne, 
Parkville, Victoria 3052, Australia 
and 

Department of Mathematics, The Australian National University 
Canberra, ACT 0200, Australia^ 

Abstract 

The restricted SOS model of Andrews, Baxter and Forrester has been studied. 
The finite size corrections to the eigenvalue spectra of the transfer matrix of the model 
with a more general crossing parameter have been calculated. Therefore the conformal 
weights and the central charges of the non-unitary or unitary minimal conformal field 
have been extracted from the finite size corrections. 
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1 Introduction 

The ABF restricted solid-on-solid (RSOS) model has been found by Andrews, Baxter and 
Forrester (ABF) in 1984 [0. It has been well known that the model provides realizations 
of unitary minimal conformal field theories [0, ^, |3J]. This has been further confirmed by 
studying the finite-size corrections to the ground state energy |5[]-[pH] (also see fl2fl-f22f 



for related works). Among these works, much effort has been focused on the ABF model 
corresponding to the unitary minimal conformal field theories. By contrast, the finite-size 
corrections to the transfer matrix of the ABF model corresponding to non-unitary minimal 
conformal field theories have received no attention. 
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The local height probabilities of the ABF model with a crossing parameter A = kn/h, 
where two relatively prime positive integers satisfy k < h, have been calculated in [p3[| . In 
this paper, with the same motivation, we repeat the consideration of the finite-size correction 
calculation of the ABF model with A = kir/h. In general the model will no longer be 
physical as there will be some negative face weights. Nevertheless, the non-unitary minimal 
conformal field theories J3J can be realized as the critical continuum of the ABF RSOS model 
with the crossing parameter A = kn/h. The model is therefore of independent interest for 
this feature. 

In [^0| an analytic method has been presented to find the finite-size corrections involving 



the central charges for the six-vertex model with a twisted boundary condition. The method 
has been successfully applied to the other models (see |27| for example). In these works only 



the central charges have been obtained. In fact the central charge and conformal weights 
together could appear in the finite-size corrections to the eigenvalue spectra of the transfer 



matrix. In this paper, following the calculation presented in p0 |, we find the finite-size 
corrections to the eigenvalue spectra of transfer matrices of the critical ABF RSOS model. 
From the corrections both the central charges and the conformal weights of non-unitary 
minimal conformal field theories are extracted. This generalizes the method presented in 



20|] to find the conformal weights of the ABF SOS model. 

We first review the ABF RSOS model and the Bethe ansatz solutions of transfer matrices 



briefly in section |1.1| . In section ^| we find an integral nonlinear equation and express the 



finite-size corrections in terms of the solution of the nonlinear equation. Then the effective 
central charges including the conformal weights are extracted from the finite-size corrections. 
A brief discussion is presented in the final section. 

1.1 Models and Bethe ansatz solutions 

The ABF RSOS model can be given by Baxter's SOS model, which was introduced in order 
to solve the eight-vertex model with the i?-matrix 

( a{u) d{u)\ 
b(u) c{u) 
c{u) b{u) 
\d{u) a(u) J 

where 



R(u) 



a(u) = e(X)e(u)H(u + A) , b(u) = Q{X)H{u)Q{u + A) , 

c(u) = H(\)Q(u)Q(u + A) , d(u) = H(X)H(u)H(u + A) . (1.2) 

The i?-matrix satisfies the Yang-Baxter equation p4|, 

R 12 (u)R 13 (u + v)R 23 (v) = R 23 (v)R n (u + v)R 12 (u) . (1.3) 

Baxter has shown in [^] that the eight- vertex model can be transferred into the SOS model, 
which is defined by the following face weights 

h(u + A) 



w{£±i,e±2,e±i,e\u) 



h(X) 



where the height £ G Z and £ is an independent parameter. The crossing parameter is A 
and the spectral parameter is u. The function h(u) is given by 

h(u) = e(0)H(u)Q(u) . (1.5) 

The face weights satisfy the following Yang-Baxter equation 

W(a, b, g, f\u)W(f, g, d, e\v)W(g, b, c, d\v-u) 

g 

= H W(f, a, g, e\v-u)W(a, b, c,g\v)W(g, c, d, e\u) (1.6) 
g 

for any integers a, b, c, d, e, /. Therefore the SOS model is an integrable system. Suppose 
that I and m are allowed spin configurations of two consecutive rows of an N (even) column 
lattice with periodic boundary conditions In+i = m N+i = m i- The elements of the row- 
to-row transfer matrix T of the SOS model are defined by 

N 

(l\T{u)\m) = J] W{lj,l j+l ,m j+u mj\u) . (1.7) 

3=1 



We recall the eigenvalues of the transfer matrix T given in [^5|] (also see |26[ for algebraic 
Bethe ansatz), 



T{u) = e isX h N (u + \\) q{U t , A) + e- isX h N (u - ±A)^±^ 
2 g(«) z g(-u) 



;i.s) 



where g(w) is defined by 



N/2 

q(u) = l[h(u-u j ). (1.9) 

These parameters «i, w 2 , • ■ • , «tv/2 are determined by the Bethe ansatz equations, 

p(« i ) = -l, j=l,2,.--,N/2 (1.10) 
where the function is given by 

^ V ; h N (u+ \X)q{u - A) 1 ; 

The ABF RSOS model is specialized by setting 

A = ktr/h and f = (1.12) 

where k and ft, are relatively prime integers (h > k > 0) and s = 1, 2, • • • , h — 1. With this 
condition (|1.12|) the face weights still satisfy the Yang-Baxter equation. The row-to-row 
transfer matrix T(u) forms the commuting family 

[T( U ), T(v)] = 0. (1.13) 



Therefore the model is integrable. The Bethe ansatz solutions ( p..8|) and (|1.10Q with the 
restriction ( 1.12| ) are the eigenvalues and the Bethe ansatz equations of the transfer matrix 
of the RSOS model 0, §. 



2 Finite-size corrections 



We consider the corresponding critical ABF RSOS model, which can be obtained by taking 
the zero elliptic nome p = 0. The elliptic function h(u) reduces to the trigonometric function 

h(u) = sin(u) (2.1) 



if p — > 0. The eigenvalues ( |1.8[ ) and the Bethe ansatz equations ( 1.10|) are still correct for 
the critical RSOS model if the function h(u) is replaced with (pT|). 

Let us introduce the new spectral variable v = iu. It is very helpful to notice that the 
eigenvalue spectra ( |1.8|) and the Bethe ansatz equations ( |1.10| ) are the same as those of 
the transfer matrix of the six vertex model with a twisted boundary condition P0| . They 



therefore can be treated similarly. The functions have to be restricted in some analyticity 
domain since all functions are i7r-periodic. It has been shown in [2(| that the following 
functions are analytic and non-zero (anz) 

h(v) ANZ in < 3m (v) < n 
q(v) ANZ in —7i < 3m (v) < 

T(v) ANZ in -A/2 < 3m (v) < A/2 , (2.2) 
and, respectively, the functions q and p satisfy 

q{v) = q(v) and p{y) = l/p(v) . (2.3) 



Because of i7r-periodic functions of face weights we can take k < h/2. Note that ( |2.2p has 
restricted the model to stay on the critical line of regime m/iv . 



2.1 Nonlinear integral equation 



Following PD|, let us introduce new functions 



MX) 



l/p{x-i\/2) = 
A(x) := 1 + a(x 



TlX 

tanh — 
2A 



N 



a[x) 



(2.4) 



The variable x may be regarded as real.[] The method presented in |2(J is to derive a set of 
relations about functions a and q and these relations lead to a nonlinear integral equation, in 
turn, the nonlinear integral equation ensures that the finite-size corrections to the eigenvalue 
spectra of the transfer matrix can be solved through dilogarithmic functions. 

The involved functions are ANZ in the strips (|2.2|) and are exponentials in asymptotic 
behaviour. The second logarithmic derivatives of the functions can be Fourier transformed, 

1 

2^ 



/(*) 



[lnf{x)\ e- lkx dx 



[In /(X)]' 



f(k) e lkx dk 



(2.5) 



4 Sometimes it is convenient to work with values of x in the upper half plane close to the real axis for 
avoiding singularities which might otherwise occur. 



where the integration path in x-plane has to lie in the analyticity strip and the real part of 
the variable of integration goes from — oo to oo. By Cauchy's theorem all other details of 
the path are irrelevant for f(k). 

We now derive a set of relations about functions a and q. Applying the Fourier transform 
to the definition (12.4) of a(x) 



a\x) 



,21sA 



, irx 
COth 2A 



A' 



h N (x)q(x - 3iA/2) 



h N (x - iA + m)q(x + iA/2 - in) ^ 

where all arguments of the functions q and h have been reduced to the analyticity strips 
( |2.2|) because of the 7ri-periodicity, then yields 



a(k) 



Nk 



+ 



Nk(l - e^ x -^ k ) 



+ ea 



T AA: 



q(k). 



1 + er Xk 1 - e _7rfc 

To solve a and q we need another relation, which is introduced by an auxiliary function 

1 + p(v 



(2.7) 



K{v) :-- 



q(v) 



(2.8) 



It is ANZ in the strip —A/2 < 3m (v) < A/2. To apply Cauchy's theorem to the Fourier 
transform of the second logarithmic derivative of h a we rewrite h a (v) in the following two 
different forms such that the arguments of q stands in the analyticity strip (|2.2|) 



K{x - iA/2) 

h a (x + iA/2) 
Applying the Fourier transform to 

e Xk/2 K{k) = 
e- Xk / 2 K(k) = 
Then they are equated yielding 

q(k) (e^ +x ^ k - e Xk ' 2 



coth 



TTX 



1 N 



A(x) 



2A 
A{x 



q(x — iA/2) a(x) 



q(x + iA/2 - ivr) 
, it follows that 
Nk 



(2.9) 



e Xk/2 q(k) + A{k) - a(k) 



1 + e- Xk 
A(k) - e (n - x/2]k q(k) . 



Nk 



-Xk 



+ a(k) + e Xk A(k)-A(k) 



(2.10) 



The equations ( |2.7| ) and ( p.!0[ ) together determine the functions a(k) and q(k), 

smh(^iik — Xk) 



a(k) 
q(k) 



2 cosh(|A/c) sinh |(7r/c — Xk) 

Nke~ nk/2 
4sinh(|7rA;) cosh(|AA;) 

p -(7T + A)fc/2 



(A(k) - e {x - e)k A{k)) 



(A(k) - e Xk A(k) 



(2.11) 



4 cosh(|A/c) sinh ^(nk — Xk) 

where an infinitesimal positive e has been introduced for the imaginary part of the argument 
of x. Transforming back to the variable x 

[lna(x)]" = [°° (K(y)[lnA)"{x-y)- K(y + ie-iX)[ln~J\"{x-y))dy (2.12) 



1 foo sinh(±7r - X)k 



where the kernel function 

K(x) := — / ; i'/" ttt e lkx dk (2.13) 

v ; 2W-oo2cosh(iAfc)sinhi(7r-A)fc V ; 

satisfies 

K(x) = K{-x) , K{x) = K{-x) . (2.14) 



The equation ( |2.12| ) is derived based on the essential anz property of the Bethe ansatz 



solution ( |1.8|) . Low- lying excitations have the same bulk behavior as the ground state. The 
only difference has been shown in || to lie in the fact that the eigenvalue functions now 
possess a finite number of zeros in the analyticity strip, which were free of zeros in the 
ground state. However, it is always possible to take an anz area in the analyticity strip 
where Cauchy's theorem can be applied [pLT]| . Therefore the equation (|2.12|) still works for 
the excited states if we change the integration path in the anz area. Integrating (|2.12|) 
twice we obtain a nonlinear integral equation 

lna(x) = I (K(y)\n.A{x — y) — K(y + \t — \\)\\iA(x — y)\dy 

+C + Dx (2.15) 

where the integral constant D = because all terms remain finite for x — > oo and another 
integral constant C heavily dependent on the branch choice of lna(x), 



C = lna(oo)-/ K(y)dy (hi A(oo) -\nA(oo] 

J —oo 



ln(cu 2 e 2isX ) - (ln(l + cu 2 e 21sX ) - ln(l + u~ 2 e 

ix — A ^ 



(2.16) 



n- A 

where the phase factor has been introduced by 

(p = sX-ilnu , u 2 = l. (2.17) 

Here we have taken a more general choice of branches for lna(x), or lna(oo) = 21n(a;e lsA ). 
The case u — 1 has been studied in f2"0[| , which corresponds to the ground state. For the 
excited states we have chosen the other branches with u ^ 1 or take 

u = e i{t - s)w (2.18) 

with the integers s,t. From the definition (|2.4j ) it follows that a goes to a e _41sA under the 
change s — >• h — s. So the same symmetry should be imposed on the equation (|2.15 ), or the 



phase factor must go to —0 under this change. It follows that the phase factor will go 
to —0 if changing s — > h — s and t — > h — k — t. Similar to the exponent s, suppose that t is 
positive. According to s = 1, 2, • • • , h — 1 we therefore take t = 1, 2, ■ • • , h — k — 1. Recalling 
the definition (|2.4|) we arrived at the nonlinear integral equation for a 



1 / \ AT! 1 nX 11T( - 

maix) = N In tanh — - + 



2A 7r — A 



+ r (K(y) lnA(x -y)- K(y + ie - iA) ln^4(a; - y)) dy . (2.19) 



This equation is exact for all finite system size and for both the ground state and the excited 
states. 



2.2 Scaling limits 

To obtain the finite-size corrections to the eigenvalue spectra of the transfer matrix we 
observe the following scaling behaviour 



lim tanh 



N 



exp (— 2e x 



(2.20) 



in thermodynamic limit iV — > oo. The function a scale similarly 

A 



a±{x) : 



A+(x) :-- 



lim a(±— (x + \nN)) , h±(x) := lna±(x) 

lim A(±-(x + \nN)) = 1 + a±(x) , lA±(x) := \nA±(x) 



(2.21) 
(2.22) 



In the scaling limit regimes the nonlinear integral equation (|2.19|) becomes 

/oo 
Ki(x - y)lA±(y)du 
-oo 

r°° m<j) 

- K 2 {x-y)lA±{y)du + 

J — oo 

poo 

h±(x) = -2e~ x + 

/oo 
K 2 (x-y)lA±(y)du- 
-oo 

where Ki^x) are defined by 

K x {x) : 
K 2 (x) : 



7T — A 

K t (x - y)lA±(y)du 

m(j) 

7T — A 



(2.23) 



7T V 7T . 



-K { —x ± i(e - A) 

7T V 7T 



(2.24) 



Let us now turn to the eigenvalues T given by ( pL.8|) . Its finite-size corrections can be 
derived from 



T(x - iA/2) = &"(x - iA) ^ + iA /2 itt) _ aiA 

g(x — iA/2) 



^2.25) 



Applying the Fourier transform to the ratio of the g-functions and taking ([2.1 1|) into account 
we have 



In 



q(x + iA/2 - m) 
q(x - iA/2) 
i 



—N 



00 sinh \(n — \)k 



e lkx dy 



+ 



2A 



oo 2k sinh \ (71k) cosh -= (Xk) 

00 In A{x-y) 



h\A(x — y) , i 
;dy + 



sinhy(y — ie) 2A J-00 sinh |(y + ie 



-dy + f c 



(2.26) 



where f c is an integration constant. Therefore the finite-size corrections to the eigenvalue 
can be expressed as 

r°° sinh ~(7r — A)A; sinhfxA;) 

lnTx-iA/2 = f c + N\nh(x-i\)-N 2 / i , rr , i „ <tt 

v 7 7 17 v 7 7-oo 2fc sinh i(7rfc) cosh l(AJfe) 

i Z" 00 ln>%) / 1 \ . „. 

+ T/ ^TTZ: - ..^ y + ilTr • ( 2 -27) 



A J-oo sin — y + ie) \N 



The scaling limit of the corrections can be done by splitting the integral into two parts, 
then inserting the variable of integration y by ±-(y + In AT) and using the scaling functions 
( p. 22 ), we obtain 



lnr(ac - iA/2) 



2 



DC 



-Nf(x - iA/2) - —ex* J lA^e'Uy 
2i r°° , / 1 

+ 



= -i\r/(x-iA/2)-^sinh™(^ f^e IA ± (v)e-«d V ^ + (i) (2.28) 

where the bulk behavior is entirely expressed by the first term and second term is the finite- 
size corrections. The integration constant f c is chosen so that f(x — iA/2) is exactly the 
bulk energy, which can be derived from the inversion relation of the face weights |7], |23f . 
Here we are only interested in the finite-size correction terms which include the conformal 
spectra. 

2.3 Conformal spectra 

The conformal spectra can be extracted from the finite-size corrections of the transfer ma- 
trix. The integral in the finite-size correction term in ( 2.28|) can be calculated by considering 
the expression 



[k±(x)] lA±(x) — 1a±(x)[lA±(x)]') dx 
+ r ([k±{x)}'M±{x) -h ± (x)[TA±(x)}') dx 

J —oo ^ ' 

2 r e~ x (lA±(x) + [lA±{x)}') dx 

J —oo ' ' 

+2 f°° e~ x (M±(x) + \jA±(x)}') dx 

J —oo ^ ' 

— ^ f°° (PM*)]' - dx . (2.29) 

71 — A J-oo v 7 



The right hand side is derived by using the nonlinear integral equation (|2.23| ). The left hand 



side can be calculated after changing the variable x to a and a and using the dilogarithmic 
function 



L(x) 



ln(l +y) \ny 
V 1+2/ 



dy . (2.30) 



Then we are able to derive 

24 ' x 



— / e~ x ^e 1A±(x)dx 

71 J-OD 



^(l(A-) + l(^V-)-M) (2.31) 

-2 



where the asymptotics of a±(oo) = (uje^ sX ^j , a±(oo) = ^e^ A ^ anda±(— oo) = a ± (— oo) 
have been read off from (|2.23|) . Finally using the well known identity 



L(z) + L(l/z) = ^ (2.32) 



the finite-size corrections in (|2.28 ) are given by the explicit expression 



lnT(x-iA/2) = -Nf{x -i\/2)- ^(c-24A)sinh^ + </— ) (2.33) 



6N\ / X \N 

or changing the variable x to v = \u = x — iA/2 

\nT(v) = -Nf(v) + ^(c-24Ay sh^ + o^ , (2.34) 
where the central charge is 



c=l- (2.35) 
tx[tx — X) 



and the conformal weights are 



Att(tt - A) 



(2.36) 



For the ground state, s = t = 1 yields A = 0. The choice ofl<s</i — 1 and 
1 < t < h — k — 1 gives the excited states. Remarkably, inserting A given by ( |1.12j ) into the 
conformal spectra we have the central charges and the conformal weights of the primary 
fields for Virasoro minimal models 

Qk 2 , . \ht - (h - k)s} 2 - k 2 

c=1 ~T7l m and A = [ tttt — ~T\ 2.37 

h(h-k) Ah(h-k) 

k < h; s = 1, 2, • • • , h — 1; t = 1, 2, • • • , h - k - 1 
for k > 1. The unitary minimal models are given by taking k = 1. 



3 Discussion 

In this paper we have obtained the conformal spectra of the non-unitary minimal conformal 
field theories from the finite-size corrections to the eigenvalue spectra of the transfer matrix 
of the critical ABF model on the regime w/n critical line with the crossing parameter ( |1.12|) . 



The method given in [T2(J] is only for calculating the central charges for the six-vertex model 
with a twisted boundary condition. In this paper it has been generalized to calculate 



both the central charges and the conformal weights. Other methods, for example, the 
thermodynamic Bethe ansatz (TBA) analysis (see ||, J28|-|35[]), exist for calculating the 



conformal spectra. The TBA relies heavily on the string hypothesis, while our method 
crucially depends on the anz property instead. However, it is an interesting problem to 
generalize the TBA method for calculating the conformal weights of the ABF SOS model. 

There is another method for calculating the finite-size corrections of transfer matrices. 
This has been shown by solving the fusion hierarchies of the ABF model. Unfortunately it 



is only for k — 1 [|Tl|] (also see |21|). It has not yet known how to find finite-size corrections 



of the transfer matrix of the ABF model for k > 1. 
Note Added 

After this work was submitted I was informed by Murray Batchelor about reference f36| 



where the authors generalize the method of |2(J to calculate the conformal spectrum of 
the six-vertex model with twisted boundary conditions. However, the finite-size corrections 
to the transfer matrix of the ABF SOS model with the more general crossing parameter 
(1.12), which is an important class of integrable lattice models corresponding to realizations 
of non-unitarity minimal conformal field theories, was not considered there. I am grateful 



to Murray Batchelor for drawing my attention to |36|. 
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